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Abstract 

We investigate the influence of the external tidal field of a dark matter halo on the dynamical evolution of star clusters 
using direct iV-body simulations, where we assume that the halo is described by a Navarro, Frenk & White mass profile 
which has an inner density cusp. We assess how varying the mass and concentration of the halo affects the rate at which 
the star cluster loses mass and we find that increasing halo mass and concentration drives enhanced mass loss rates and 
in principle shorter cluster disruption timescales. In addition, we examine disruption timescales in a three-component 
model of a galaxy (bulge, disk and dark matter halo) and find good agreement with results based on an empirical model 
of the Galactic potential if we assume a halo mass of ~ 10 12 M Q . In general, dark matter halos are expected to contribute 
significantly to the masses of galaxies and should not be ignored when modelling the evolution of star clusters. We extend 
pur results to discuss how this can have a potentially profound effect on the disruption timescales of globular clusters, 
suggesting that we may underestimate the rate at which primordial globular clusters are disrupted. 

'Key words: globular clusters: general, galaxies: halos, dark matter, methods: n-body simulations 
PACS: 98.10,+z, 98.35.Gi 
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■1. Introduction 

i Globular clusters are compact stellar systems, each con- 
taining of order 10 6 stars, that orbit around the centres 
of galaxies out to large radii. Our Galaxy hosts around 
200 globular clusters and observations show that almost 
all galaxies host these systems, with giant ellipticals hav 



ing t he largest (relative) population (jBrodie and Strade 
200(J ). Stellar population studies have reve aled that glob 



ular clusters have ages up to ~ 13 Gyrs (e.g. Hansen et all 
l2002tlChabover and Kraussl , l2002[h This implies that they 
formed within ~ 1 billion years of the Big Bang and so they 
represent fossil records of the earliest epoch of galaxy for- 
mation, which has led to much interes t in their potential as 
probe s of the high redshift Universe (jBrodie and Straderl . 



However, we have as yet no compelling theory for glob- 
ular cluster formation and evolution within a cosmological 
framework and this has important consequences for the 
kind of questions we can address using globular clusters 
and the strength of conclusions that we can draw. One 
particularly important consequence concerns the disrup- 
tion of primordial globular clusters and the relation of the 
high redshift globular cluster population to the popula- 
tion that we observe at the present day. All star clusters 
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lose mass over time and this depends on a number of fac- 
tors including initial conditions (e.g. initial mass profile 
of the star cluster, initial mass function of the stars, ini- 
tial binary fraction), internal processes within the cluster 
(e.g. mass segregation, stellar mass loss), the orbital pa- 
rameters of the cluster and the form and nature of the 
external tidal field within which the cluster orbits. The 
globular cluster population that we observe around galax- 
ies at the present day has survived a Hubble time, but 
it is interesting to ask whether or not this population 
is representative of the primordial globular cluster pop- 
ulation at high redshift. If not, this would imply that 
many primordial globular clusters disrupted on relatively 
short timescales, and that the clusters that survived to the 
present day arc in a sense atypical. This introduces cer- 
tain caveats where the present day population is used to 
probe, say, the effici e ncy o f star formation at high rcdshifts 
( Spitler and Forbes . 2008 ) and the epoch of cosmological 
reionisation (jMoore et al. , 2006^ 1 . 

In the standard picture of galaxy formation, galaxies 
are embedde d in massive v i rialise d structures or halos of 
dark matter (jSpringel et al. . 120061) . These dark matter ha- 



los assemble over time hierarchically, continuously growing 
via accretion of dark matter and merging with other ha- 
los from high redshift to the present day. The subset of 
globular clusters that have survived to the present day 
has evolved in a time-dependent and at times violently 
changing potential provided by its host galaxy and dark 
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matter halo. Previou s studies that have investigated star 
clust e r mass loss ( e.g. Giersz and Heggiel 19971 : IVesperini , 
1998t iBaumgardtl . l200lt iHurlev et all booU T^ with mod- 
els ranging in mass from open clusters to small globular 
clusters - have assumed simplistic models for the Galaxy's 
tidal field by treating the Galaxy as a central compact 
source or by using Oort constants that are empirically 
meas ured from stellar motions in th e Solar neighbourhood 
(e.g. Binnev and Tremainei . Il987 ) . Mashchenko and Si lis 
( 2f)05h looked at the early evolution of globular clusters 
formed in dark matter halos but employed softened gravi- 
tational potentials for both the stars and dark matter. As 
such, the explicit presence^] of a dark matter halo has been 
neglected to date in direct iV-body models of star cluster 
evolution. 

Our goal is to understand how the presence of a dark 
matter halo affects the internal characteristics of star clus- 
ters in general - and globular clusters in particular - and 
their disruption rates. For the former our interest is in 
signatures such as the ratio of the core to half-mass radii 
and the velocity dispersion, which may be compared to 
observations of star clusters. However, in this paper, the 
first in a series, we start by focusing on the rate of mass- 
loss from model star clusters and what this can tell us 
about the expected lifetimes and disruption of globular 
clusters. Our approach is to use a direct iV-body code that 
treats the internal processes of a star cluster in as detailed 
a manner as possible. We include a spectrum of stellar 
masses and stellar evolution from the outset as the effect of 
t hes c on star cluster evolut ion is relatively well understood 
(e.g. de La Fuente Marcos! . 1996t Baumgardt and Makinol . 
2003rlHurlev et all 120041 ). In our treatment of the external 
potential and the orbits of the model clusters within this, 
we take a more cautious approach. We start with the sim- 
plest case, circular orbits in a variety of static potentials, 
as this facilitates comparison with previous iV-body stud- 
ies and provides the foundation for our study. In the future 
we will expand the models to include eccentric orbits and 
time-varying potentials. By necessity we also start with 
models more comparable in size to open clusters than glob- 
ular clusters but this will also increase as we progress. In 
this way we build a consistent and realistic picture of star 
cluster evolution and survival. 

We assess the influence of halo structure on mass loss 
rate using the functional form for th e mass profi le pro- 
posed by iNavarro. Frenk. and White! (|l996l . Il997t here- 
after NFW). The NFW profile has been found to provide 
a good description of the spherically averaged dark matter 
distribution within halos in dynamical equilibrium in cos- 
mological iV-body simulations. The characteristic feature 
of the NFW profile is that it has a central density cusp, 
p oc 



so that the central density is divergent. We also 
construct a three-component model of our Galaxy consist- 



x We use explicit because models of the Galactic potential based 
on Oort constants implicitly include the contribution of the Galaxy's 
dark matter halo. 



ing of a bulge, disk and NFW halo and we compare the 
rate at which clusters lose mass in this model against the 
rate measured in the fiducial model of the Galactic tidal 
field, based on Oort constants. 

The format of the paper is as follows. In the next sec- 
tion we describe the NFW profile in more detail, while 
in section [3] we describe the model we have adopted for 
our star cluster. In section 3] we present the main results 
of our study, first considering the evolution of an isolated 
cluster; then the evolution of the cluster in external galac- 
tic potentials; then evolution in the potential of a NFW 
halo; and finally evolution in our three-component model 
of a galactic potential, which consists of a bulge, disk and 
NFW halo. Finally we discuss the main findings and im- 
plications of our study in section [5l 

2. The NFW Mass Profile 

Dark matter halos that form in cosmological simula- 
tions arc relati vely complex structures - the y are generally 
aspheri cal (e.g. Bailin and Steinmet3 . 12005! ) and asymmet- 
ric (c.g. lGao and Whit el. l200ll , they have no simple bound- 
ary (e.g. Prada et all 12006^ and they contain a wealth of 
small scale structure ( Gao et al 1 12004 . Despite this com- 
plexity, however, it is conventional to identify a halo as a 
spherical overdense region, typically of order 100 times the 
mean density of the Universe. The mass enclosed within 
this spherical overdensity defines the virial mass of the 
halo. 
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A vir Pcrit ^Vir " 



(1) 



Here p cr ; t — 3H 2 /8irG is the critical density of the Uni- 
verse and r v ; r is the virial radius, which defines the ra- 
dial extent of the halo. A v ; r , the virial overdensity crite- 
rion, corresponds to the mean overdensity at the time of 
virialisation in the spherical collapse m odel, the s i mples t 



analytic model of halo formation (cf. lEke et a~ . 1996h . 



Depending on redshift and cosmological parameters, A v j r 
varies between ~ 100 and ~ 200. 

The mass profiles of dark matter halos in dynamical 
equilibrium forming in cosmological A-body simulations 
can be relatively well described by the NFW profile, 



p(r) 



Pcrit^c 



r/r s (1 + r/r s ) 2 



(2) 



where p cr it is the critical density of the Universe, r s is the 
scale radius and i5 c is the characteristic overdensity. The 
characteristic overdensity 5 C is itself a function of r s , 



S c = 



Avir 



where c 



3 ln(l + c) - c/(l + c) ' 
/r s is the concentration. 



(3) 



Because the scale radius r s and the characteristic over- 
density S c arc related, cquation.r5]can be rewritten in terms 
of a single parameter, the concentration c = r v - lr /r s , such 



2 



Model Description 

1 Isolated cluster 

2 Standard Galactic tidal field 

3 Two-component Galaxy (disk and bulge) 

4 NFW halo 

5 Full Galaxy (disk, bulge and NFW halo) 

Tabic 1: Overview of the models used in the simulations in terms of 
the treatment of the external tidal field. 



that for a fixed concentration, the local density depends 
only on the normalised radius r/r V j r . Cosmological simula- 
tions have shown that virial mass M v - a and concentration c 
are correlated, such that t he concentration increases as the 



virial mass decre a ses (e .g iBullock et all 120011 : lEke et al 



2007). As a rule of thumb, a typical 



20011: iNeto et al. 

M vir =10 10 (10 12 , 10 14 ) M Q halo has a concentration c ~ 15 
(10,5). 

The mass enclosed within a radius r is obtained from 
equation [2J 

M(< r) = ATTp cHt 5 c r^g{r) 

where 

r/r s 



(4) 



g(r) = ln(l + r/r s ) 



1 + r/r s 



(5) 



This results in a radial acceleration directed towards the 
centre of the halo of 



a(r) = -AwGp C n t S c r s g(r) (r/r 8 )~ 



(6) 



It follows that the gradient in the radial acceleration (and 
consequently the force) is given by 



da 
Or 



-^- = -2 



a{r) 



AnGp(r) 



(7) 



This provides a measure of the tidal force across the cluster 
and in the absence of an extended dark matter halo the 
tidal force would reduce to the familiar result for a point 

mass, 



da 
dr 



— = -2 



GM(r) 



(8) 



3. The cluster model 



To study the evolution of a star cluste r we used the 
iV-body code NBODY6 (see lAarsethl . 120031 for a full de- 
scription). It uses Hermite integration with individual 
time-steps and does not use softening in the force equa- 
tion. Stella r evolution is includ ed using the alg orithms 
provid ed by iHurlev et al. ( 2000h as described in Hurley 
( 2008h . The code also allows for binary formation and 
evolution. There are also options to include an external 
Galactic potential with the simplest option being a stan- 
dard Galactic potential at 8.5 kpc based on the local Oort 



constants (hereafter the standard Galactic tidal field). In 
most resp ects NBODY6 i s very similar to it s sister code 
NBODY4 (jAarsethl . Il999l ; IHurlev et ail . l200lh except that 
the la tter interfaces with special-purpose GRAPE hard- 
ware ([Making et all 120031 ) to speed up the calculation of 
the gravitational forces between the stars. NBODY6 op- 
erates on standard hardware and our simulations are per- 
formed on the Swinburne supercomputeiQ. 

Within NBODY6 there is the possibility of building a 
multi -component galaxy consistin g of bulge (point mass), 
disk ( Miyam oto and. Nagail . Il975h and logarithmic poten- 
tial ( AarsethT 2003^ 1 . To this we have added the possi- 
bility of replacing the logarithmic potential with a po- 
tential based on the NFW-profile. This three-component 
galactic potential thus consists of a bulge, disk and NFW 
halo. In all other respects the NBODY6 software is the 
same as the version publicly available for download at 



http : //www. ast . cam. ac .uk~sverre/web/pages/nbody .htm 

with input options and parameters set as for the defaults 
sugges ted in the included manual and described bv lAarsethl 
( 2003h . 

All our models start with si ngle stars o nly a nd these 
are distributed according to the IPlummerl (| 1 9 1 lh density 
model, an n = 5 polytrope. We assume the stars are in 
virial equilibrium when setting the initial positions and ve- 
locities and take the initial virial radius to be 3 pc which 
sets the length scale of the simulation. It should also be 
noted that the Plummer profile formally extends to infi- 
nite radius so a cut-off at 10 times the half-mass radius is 



appli ed to deal with rare cases of large distance (jAarsethl . 
20031) . To save computing time we did not run full globular 
cluster models (i.e. TV = 10 5 — 10 6 ). Instead we focus on 
models of 1000 stars. For one galaxy scenario we also ran 
N = 16 000 simulations in order to investigate the scaling 
of the results with N. The initial masses of the stars are 
distributed according to a Salpeter initial mass function 
(IMF: n(M) oc M -2 3 ) with a minimum and maximum 
stellar mass m m i n = 0.3 M Q and m max = 30 M Q respec- 
tively. The average stellar mass is 1 M Q . In general our 
clusters are placed on circular orbits in the stellar disk (if 
present) at a radius of 8.5 kpc from the Galactic centre. 
This facilitates comparison with previous studies. How- 
ever, we do briefly explore the effect of orbital distance 
by placing some clusters on orbits at radii of 4.5, 17 and 
34 kpc. An overview of our models used in terms of the 
treatment of the external tidal field is shown in Table [TJ 

4. Results 

4-1- Star Cluster Evolution in Isolation 

Our first model (Model 1) is a cluster in isolation. Its 
dynamical evolution depends only on internal processes 
such as mass segregation and stellar evolution. The evolu- 
tion of the total mass of such a cluster is shown in Figure 



2 scc http:/ /astronomy.swin. edu.au/supercomputing/ 
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Figure 1: The evolution of the total mass of the cluster in isola- 
tion (Model 1, solid), in the standard Galactic tidal field (Model 2, 
dashed) and in the two-component galaxy (Model 3, dotted). All 
models start with N = 1000. 



[T] by the heavy solid curve. In NBODY6 the cluster size 
is limited by a tidal radius or cut-off radius. For the case 
of an isolated cluster it is not meaningful to define a tidal 
radius and so we adopt a cut-off radius of 10 times the 
cluster scale length r s , 



1 gMgrt 

2 $ : 



(9) 



with Mtot the total cluster mass and $ the potential en- 
ergy. Note that r s is distinct from the NFW scale radius 
discussed in section Stars with cluster-centric radii in 
excess of twice the cut-off radius are considered to have 
left the cluster and contribute to the mass loss. The ini- 
tial cut-off radius is 23 pc, but this increases as the cluster 
expands to about 100 pc after 3Gyr. 

Mass segregation causes the system to lose stars early 
on. Less massive stars gain kinetic energy through equipar- 
tition during encounters with more massive stars in the 
core and are more likely to escape. Stellar evolution leads 
to the loss of even more mass in the early stages. Massive 
stars evolve more rapidly than low mass stars and hence 
lose mass more quickly in the form of stellar winds. Neu- 
tron stars that form in a supernova will get a velocity kick 
and generally disappear from the system. 

4-2. The Standard Galactic Tidal Field Model 

Model 2 introduces the standard Galactic tidal field 
(SGTF) based on the local Oort constants (A = 14.4 km/s, 
B = —12.0 km/s). These are derived empirically from 
the observed motions of star s in the Solar neighbourhood 
(jBinnev and Tremaind . 119871 ). As such they naturally ac- 
count for the local effect of the total mass of the Galaxy 
(including its dark matter halo) and they strictly apply to 
orbits within the region of the Solar neighbourhood. From 
the local circular velocity we can infer that the enclosed 



mass at an orbit of 8.5 kpc (the distance from the Sun to 
the Galactic centre) is 9 x 1O 1O M0. As the dashed curve 
in Figure [1] shows this has a profound effect on the mass 
loss of the cluster. There is an almost continuous mass loss 
leading to rapid disruption of the cluster: simulations are 
halted when the cluster has less than 25 stars left. Here 



we define the tidal radius as the iKingl (jl962[ ) tidal radius, 

q 1 MC 9 , > 

* 3M G G v ' 

with Mc and Mq the masses of the cluster and the galaxy 
respectively and Rq the distance between the galactic cen- 
tre and the centre of the cluster. We set Rq = 8.5 kpc (the 
distance from our Sun to the Galactic centre). This gives 
a tidal radius of ~ 16 pc, smaller than the cut-off radius 
used in Model 1. This has a bearing on what we com- 
pute as the mass loss rate, which will be lower in the case 
of the isolated cluster where the cut-off radius is higher. 
However, stellar evolution effects dominate in the isolated 
cluster so the general comparison holds. 

4-3. The Two-Component Model - Bulge & Disk 

Model 3 introduces the two-component model of the 
Galaxy, consisting of a bulge a nd disk. The disk is mod- 
elled as a Miyamoto-Nagai disk ( Miyamoto and Nagail . [l975l ) 
and the bulge is a point mass at the centre of the galaxy. 
We use a mass of 5 x 10 10 M© for the disk, which has a 
radius of 40 kpc, and 1.5 x 10 10 M for the bulge, as sug- 
gested bv lXue et ail (I2008T ) for the Milky Way. 

We adopt a tidal radius of 16 pc, equal to the King ra- 
dius (sec Equation [TU1) of Model 2. This is kept fixed for 
all subsequent models at an orbital distance of 8.5 kpc, in- 
cluding the models where we add a dark matter halo. Thus 
we are effect i vely e mploying a tidal cut-off as described by 



Trenti et all (|2007f ) 



The dotted line in Figure [T] shows the evolution of the 
cluster mass with Model 3 acting as the external tidal 
field. Model 3 does not produce a tidal field as strong as 
the SGTF of Model 2, but this is to be expected because 
the mass of the combined bulge and disk within 8.5 kpc is 
smaller than the enclosed mass one infers for the SGTF 
of Model 2, which is 9 x 10 10 M Q . Recall that Model 2 is 
based on an empirical estimate of the Galactic tidal field 
and so implicitly includes the contribution of the Galaxy's 
dark matter halo, whereas Model 3 has no halo contribu- 
tion. 

The two-component model is attractive because we can 
fully determine how the Galactic potential is modelled and 
explore different sets of assumptions, but it neglects the 
presence of dark matter halo component to the potential 
which we expect to be important. In contrast, the SGTF 
model provides a realistic model of the Galactic poten- 
tial in the Solar neighbourhood, being derived from em- 
pirical measurement, but the conditions under which it is 
applicable are restrictive. If the two-component model is 
to be made more realistic, then it must become a three- 
component model (a bulge, disk and halo), in which case 
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Figure 2: The total mass evolution of a cluster of 1000 stars in 
an NFW halo at an orbital distance of 8.5 kpc from the halo cen- 
tre (Model 4). We show here three different virial masses: 10 9 Mq 
(green), 1.5 X 10 12 M (blue) and 10 14 M Q (red); each with three 
different concentrations: c = 5 (solid), 10 (dashed) and 15 (dotted). 
The solid black line is Model 1, the cluster in isolation, and is drawn 
here as a reference. 



we expect the mass loss rates for Models 2 and 3 to con- 
verge. In the next subsection we look at the influence of 
a single component model - that of a NFW dark matter 
halo on the mass loss rate - before introducing the three- 
component model in the final subsection. 

44. The NFW Halo 

We now consider a set of models in which our star clus- 
ters orbit in a NFW dark matter halo (Model 4). For now 
we ignore the stellar component of the galaxy because we 
want to investigate the influence of changing the NFW 
parameters of mass and concentration. We use a model 
grid of 6 different halo virial masses (10 9 M Q , 1O 1O M0, 
10 n M o , 1.5xlO 12 M , 10 13 M Q and 10 14 M Q ) and 3 con- 
centrations (5, 10 and 15), thus 18 models in total. Here 
we define the virial mass as in equation [T] with an overden- 
sity criterion of A v i r = 200. Each cluster follows a circular 
orbit at a distance of 8.5 kpc from the centre of the halo. 
This is achieved by giving the cluster an initial circular 
velocity equal to the halo velocity at that distance, 



V c = s/(GM(< r)/r). 



(11) 



Our main set of simulations start with 1000 stars with a 
Salpctcr IMF in the mass range 0.3 — 30.0 Mq and total 
mass of 1000 M Q . Figure [2] shows the evolution of the mass 
of the cluster. For clarity we only show 9 of the models rep- 
resenting the lowest mass, intermediate mass and highest 
mass halos. The mass loss is initially the same for all mod- 
els since in this stage it is dominated by stellar evolution 
of the most massive stars. After about 100 Myr clusters in 
the most massive halos keep losing mass at a faster pace 



Mass 
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%r vir 


V c 


dV c /dr 


N 




(M ) 






(km/s) 


(km/s/kpc) 




(Myr) 


1.0 x 10 9 


5 


52 


12.46 


0.68 


1000 


853 


1.0 x 10 9 


10 


52 


13.27 


0.55 


1000 


916 


1.0 x 10 9 


15 


52 


13.66 


0.49 


1000 


851 


1.0 x 10 10 
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24 


17.77 


1.27 


1000 


859 


1.0 x 10 10 


10 


24 


20.79 


1.17 


1000 


921 


1.0 x 10 10 


15 


24 


22.42 


1.08 


1000 


868 


1.0 x 10 11 


5 


11 


22.68 


2.00 


1000 


946 


1.0 x 10 11 


10 


11 


29.34 


2.15 


1000 


826 


1.0 x 10 11 


15 


11 


33.51 


2.15 


1000 


828 


1.5 x 10 12 


5 


4.6 


26.75 


2.75 


1000 


960 


1.5 x 10 12 


10 


4.6 


38.21 


3.52 


1000 


614 


1.5 x 10 12 


15 


4.6 


46.69 


3.94 


1000 


483 


1.5 x 10 12 


5 


4.6 


26.75 


2.75 


16000 


3749 


1.5 x 10 12 


10 


4.6 


38.21 


3.52 


16000 


2746 


1.5 x 10 12 


15 


4.6 


46.69 


3.94 


16000 


2152 


1.0 x 10 13 


5 


2.4 


28.45 


3.10 


1000 


825 


1.0 x 10 13 


10 


2.4 


42.60 


4.34 


1000 


490 


1.0 x 10 13 


15 


2.4 


54.02 


5.19 


1000 


324 


1.0 x 10 14 


5 


1.0 


29.62 


3.36 


1000 


587 


1.0 x 10 14 


10 


1.0 


45.90 


5.03 


1000 


314 


1.0 x 10 14 


15 


1.0 


59.99 


6.37 


1000 


190 



Table 2: Disruption times - defined as the time after which the 
cluster has lost half of its initial mass — for clusters orbiting in NFW 
halos of different halo mass and concentration. We also show the 
percentage of the virial radius at which the cluster is orbiting at 
8.5 kpc, the circular velocity and the gradient in the circular velocity 
at that radius, and the initial number of stars in the model. 



than the clusters in the less massive halos. In these mas- 
sive halos the effect of concentration is also clearly visible. 
Higher concentrations mean faster disruptions - more of 
the halo mass is concentrated within the cluster's orbit re- 
sulting in a greater tidal force on the cluster. The effects of 
concentration become negligible for the least massive ha- 
los, and it is also noteworthy that mass loss is very similar 
for the low- and intermediate-mass halos. 

To quantify the disruption of a cluster we have defined 
the disruption time scale as the time it takes for the cluster 
to lose half of its mass. These timescales are tabulated 
in Table [H Generally, clusters orbiting in more massive 
halos have shorter disruption times than those orbiting in 
less massive halos. In massive halos the concentration has 
a large effect on the disruption time while in the lower 
mass ranges the disruption is independent of mass and 
concentration. Variations in this regime are dominated by 
statistical effects. 

We can understand why this might be by considering 
the tidal field across a cluster and determining how large 
it needs to be for the typical speeds of stars in the clus- 
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Figure 3: The mass loss for the standard NFW model (Af v ir = 
1.5 X 10 12 Mq, c = 10) at different orbital distances: 4.5 kpc (solid), 
8.5kpc (dashed), 17kpc (dashcd-dotted) and 34kpc (dotted). All 
models start with N = 1000. 
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Figure 4: The total mass evolution for a cluster at an orbital distance 
of 8.5 kpc within a full Galactic model (Model 5). Labels are as in 
Figure [5] except the solid black line now represents Model 2, with 
the tidal field based on the Oort constants. All models start with 
N = 1000 



ter to exceed the velocity dispersion. If a is the velocity 
dispersion of the cluster and R c is its radius, then on a 
timcscalc / x R c /a the net velocity of a star across the 
cluster is approximately da/dr x 2 R c x / x R c /a. Here 
the factor / will be greater than unity. If the test particle 
is to escape the cluster, 4 x / x (V c /R) 2 x R^/a must be 
comparable to a, where V c is the halo circular velocity at 
halo-centric radius R. This reduces to the condition that 



(12) 



For the star cluster mass we have considered, a ~ lkm/s 
and Rc/R ~ 0.001, so V c ~ 500km/s if / = 1, or in other 
words in heavy halos clusters disrupt very efficiently For 
less massive halos, we require / ~ 100, and so internal 
processes rather than the external tidal field will be im- 
portant in driving the mass loss rate. 

The model with virial halo mass of 1.5 x 10 12 M Q and 
concentration c = 10 resembles closest a model for our 
own Milky Way halo and from now on we will call it the 
standard NFW model. We ran N = 16 000 simulations 
for the models with this halo mass and the three differ- 
ent concentrations. The results are shown in Table [5] for 
comparison. The tidal radius for these larger models is 
40 pc, as given by equation [TOJ We see the same general 
trend of decreasing disruption time with increasing con- 
centration as seen for the N = 1000 models. It has been 



shown previously by iBaumgardtl (|200lh that the disrup- 
tion time of a star cluster can be expected to scale by the 
initial half-mass relaxation timescale and an N dependent 
factor of (log(O.lliV) /N) 1/A (see also lTrenti etaill2007t ). 



This is based on simulations of star clusters within a stan- 
dard Galactic tidal field. The disruption timescales of our 
N = 16 000 models exceed those of their N = 1000 coun- 



terparts by a factor of 4—4.5 while the half-mass relaxation 
timescale is a factor of two greater on a verage in the larger 
models. Thus we do not agree with the IBaumgardtl ([2001) 
scaling relation by almost a factor of three. However, the 
use of both a tidal cut-off and a non-standard tidal field 
will certainly play a role and we will attempt to quan- 
tify this in future work after performing simulations for an 
expanded range of N. 

For the standard NFW model we also looked at the ef- 
fects of orbital distance. Figure [3] shows the mass evolution 
for this standard model at four different orbital distances. 
A smaller orbital radius leads to faster disruption. Note 
that the tidal radius remains fixed at 16 pc for these mod- 
els. 



4-5. The Three- Component Galaxy - Bulge, Disk & Halo 

In Model 5 we added the two-component model Galaxy 
(Model 3) to the NFW models (Model 4), producing a 
three-component model consisting of bulge, disk and dark 
matter halo. Figure 0] shows the mass evolution for this 
model. The presence of the bulge and disk results in more 
rapid mass loss and consequently a shorter disruption time 
for the cluster. This is to be expected because of the 
addition of mass. The relative effect is greatest on the 
lowest mass halos because the size of the stellar compo- 
nent of the Galaxy is larger relative to the mass of their 
halo than is the case for larger halos. Comparing this to 
Model 2 (external field based on the local Oort constants) 
we see that it can be well reproduced by representing the 
galaxy as a disk, bulge and dark matter halo with a mass 
of 1.5 x 10 12 M Q and concentration c = 15. The enclosed 
mass of the dark matter halo at an orbital distance of 
8.5 kpc is 9.4 x 10 10 M . This is of order 20% larger than 
what one infers from the Oort constants. All the models 
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so far have assumed orbits in the disk. However we have 
looked at a model whose orbital is perpendicular to the 
plane of the disk; this resulted in an enhanced disruption 
rate because of "disk shocking" : the cluster experiences a 
sudden change in the potential as it goes through the disk 
resulting in faster mass loss. 

5. Summary 

We have examined the importance the presence of a 
dark matter halo has on the rate at which a star cluster 
orbiting in its tidal field loses mass. For our halo model we 
adopted the NFW profile, which provides a good descrip- 
tion of the mass profiles of dark matter halos forming in 
cosmological simulations. The NFW profile can be char- 
acterised by a single parameter, its concentration, which 
provides a measure of how rapidly the logarithmic slope 
changes with radius and consequently how rapidly the ra- 
dial force changes with radius. As expected, star clusters 
orbiting at a fixed physical radius lose mass at a higher 
rate in more massive halos and in more concentrated ha- 
los. 

We also investigated a three-component model of the 
Galactic tidal field, consisting of a stellar bulge and disk 
and a NFW halo and compared the mass loss rates of star 
clusters orbiting in this potential against mass loss rates of 
clusters orbiting in a standard Galactic tidal field (SGTF). 
The SGTF is based on empirical estimates of the Galac- 
tic tidal field in the Solar neighbourhood based on stellar 
motions, and as such provides a measure of total mass 
within the Solar radius, implicitly including our Galaxy's 
dark matter halo. We showed that star clusters orbit- 
ing in three-component models lose mass at a comparable 
rate to those orbiting in SGTF models if the halo mass is 
1.5 x 10 12 M Q and concentration c = 15. These numbers 
are reasonable and consistent with what one would expect 
for Galactic-type dark matter halos in cosmological simu- 
lations. By building multi-component models of galactic 
potentials, one can explore many questions of cluster evo- 
lution that have not necessarily been possible before (such 
as the disruption timescalcs of clusters on eccentric orbits) . 

Our results arc interesting because they highlight the 
importance of including an underlying dark matter halo 
in TV-body simulations of long-term globular cluster evo- 
lution. Previous studies have neglected to do this, but as 
we have shown, the halo can have a profound effect on the 
lifetime of the star cluster. However, in the centre of the 
galaxy baryons will dominate the mass contribution. In 
future work we will also look at the effect of the halo mass 
on the internal properties of star clusters as they evolve. 
This is an im portant step in t he spirit of the MODEST 
collaboration (jSills et all [2.003) to increase the realism of 
star cluster simulations. Of course, our calculation is ide- 
alised; we assume a static NFW halo - one that does not 
evolve over time - and our star clusters remain on a circu- 
lar orbit at a fixed radius over a Hubble time. As we have 
noted in the introduction, we would expect a realistic halo 



to have had an active and at times violent mass assembly 
history, which should be an important additional driver of 
mass loss in our clusters. Furthermore, it is unclear what 
effect the assembly of the galaxy would have on the struc- 
ture of the underlying dark matter halo. Nevertheless, our 
general conclusions hold and the effect of non-circular or- 
bits and an evolving halo will be qualified as we continue 
our study. 

Our results raise a number of interesting questions. 
The favoured cosmological model predict s that galaxies 



shoul d live in cuspy dark matter halos (jSpringel et al 
Previous studies of globular cluster evolution have 



been of clusters that evolve in non-cuspy potentials. In- 
ferences have been drawn from these simulations of gen- 
eral processes such as core-collapse and mass-segregation 
and observations of globular clusters, which the simula- 
tions seem to do a good job of reproducing. We do not 
appear to need cuspy dark matter halos to produce these 
these trends in simulated globular clusters that are con- 
sistent with observational data, so does this imply that 
cuspy dark matter halos are not needed at all? This is 
not straightforward to answer, because it is not under- 
stood how the assembly of the galaxy has affected its dark 
matter halo, and whether this has effectively wiped out the 
cusp (as would appear necessary from rotation curve stud- 
ies of dark matter dominated galaxies). However, if cusps 
are robust in the presence of the growing galaxy, then it 
may have implications for our understanding of globular 
cluster evolution. 

More speculatively, the enhanced rate of disruption of 
clusters in cuspy halos has interesting implications for their 
use as probes of galaxy formation. Understanding how 
efficiently globular clusters form should provide us with 
insight into the efficiency of galaxy formation at high red- 
shifts. The number of old globular clusters today could tell 
us how plentiful the sites of globular cluster formation were 
in the past, which in turn tells us about the efficiency with 
which gas cooled and formed molecular clouds. These sites 
may be regulated by the location (in the disks of gas-rich 
proto-galaxics, in mergers between gas-rich proto-galaxics) 
or the ambient radiation field (cosmological rcionisation). 
If globular clusters are disrupted more efficiently in cuspy 
halos, then we could be missing a large fraction of globu- 
lar clusters thus misinterpreting observations and drawing 
even more uncertain conclusions. 

These are interesting questions and ones that we aim 
to revisit in coming papers. 
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